A Hagedorn temperature in post-inflationary dynamics? 
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We consider the possibility of regarding the maximum temperature T max of a plasma of inflaton 
decay products as a Hagedorn temperature connected with the transition from an inflaton-dominated 
epoch to a radiation-dominated universe. Employing a "bottom-up" approach to inflaton decay, we 
show that the limiting temperature Tk ~ m^, as obtained from a numerical analysis of the evolution 
equations, can also be understood theoretically in a Hagedorn picture of the phase transition. 



m : 
o . 
o ■ 

(N ; 
Or 

m : 

<n : 
<n . 

> : 

in . 
i> ■ 

o\ : 
o . 

m ■ 

o ; 

^ : 

a: 
i 

<D . 
^ ■ 



X 



The concept of Hagedorn temperature in statisti- 
cal systems is commonly invoked as a scale that signals 
the onset of a phase transition 0. While it does not 
directly relate to the usual thermodynamic definition of 
temperature, it can be regarded as a limiting temper- 
ature where the population of energy states begins to 
exceed the Boltzmann weight, implying that the original 
degrees of freedom strongly overlap in momentum space 
and lose their identity. One example of this is the melting 
of hadronic and glueball states in relativistic heavy-ion 
collisions, where the Hagedorn temperature Th can be 
derived from the Bag model 3] or a stringy model [J] of 
hadrons, and can be used as an estimate of the decon- 
finement temperature, though the exact equivalence of 
the two is hard to establish, at least in four-dimensional 
gauge theories p|. However, the idea of a number den- 
sity of states that grows exponentially with temperature 
is quite natural coming from a theory with confining 
forces (like Regge theory), and may be useful when the 
exact nature of these microscopic forces is far from clear. 
A closely related example is Hagedorn-type behavior in 
some string theories at weak coupling 2] . In the context 
of early string dominated cosmology, a Hagedorn-type 
phase could be associated with open strings attached to 
D-branes which has recently been shown to provide a 
negative pressure necessary to drive inflation at high tem- 
peratures near the string scale |f|. This is accompanied 
by an exponential growth in the number of string modes 
with energy. 

In this letter, we explore the possibility to have a 
Hagedorn limit of a quite different nature in cosmology, 
namely, in post-inflationary dynamics that culminate in 
reheating of the universe. Inflation is a well known con- 
tender for making the Universe flat, besides generating 
adiabatic density perturbations with an almost flat spec- 
trum 0- However, inflation must end in order to provide 
a thermal bath of radiation in order to facilitate synthe- 
sis of light nuclei. The standard paradigm is that the 
inflaton field, 0, enters the oscillatory phase and decays 
perturbatively into radiation quanta, which subsequently 
thermalize completely through interactions, yielding a re- 
heat temperature T Te u, when the universe is completely 
radiation dominated [q. It has also been proposed that 



an initial non-perturbative decay of the inflaton, dubbed 
as preheating, is a plausible scenario [ljj. Neverthe- 
less there has to be a subsequent stage of thermaliza- 
tion after preheating quite similar to perturbative reheat- 
ing Yet another interesting possibility is to have a 
non-topological defect formation from the fragmentation 
of the inflaton condensate, such as Q-balls where 
reheating occurs via surface evaporation of Q-balls. 

Irrespective of a particular scenario the final reheat- 
ing temperature plays a very crucial role in the early 
universe, for e.g. see |l2j . The reheat temperature 
can be constrained by bounds coming from thermal 
and non-thermal gravitino production during reheat- 
ing/preheating in supersymmetric theories |l3|. Typi- 
cally, T re h < 10 9 GeV is required, for gravitino mass 
m 3/2 ~ 100 GeV, in order to avoid constraints coming 
from successful predictions of Big Bang Nucleosynthesis 
(BBN), for e.g. see Q. 

The reheat temperature is not, however, the maximum 
temperature of a thermal plasma, which may be two 
orders of magnitude higher (but less than m</,) 0, Hif . 
Physically, this maximum temperature, denoted by T max , 
arises due to the competition between early entropy pro- 
duction within light relativistic degrees of freedom at the 
earliest stages of inflaton decay and continued expansion 
of the universe which dilutes the energy density. The 
main objective of this work is to show that T max may be 
regarded as a limiting temperature, if we trace the infla- 
ton decay process backwards in time, up to a point when 
coherent inflaton production from the radiation fields be- 
gins. We take into account thermal factors that capture 
the essential feature of Bose enhancement when the in- 
flaton population builds up. We support our theoretical 
conjecture about this limiting temperature as a Hagedorn 
temperature with a numerical analysis of the evolution 
equations for matter (inflaton) and radiation energy den- 
sities run backwards in time. 

After inflation ends, the energy density stored in coher- 
ent oscillations of the inflaton field is transferred to radia- 
tion quanta over a period of time. Several features of this 
transition, such as thermalization associated heavy- 
particle production and non-perturbative pro cesses 
(preheating) have been widely investigated [lfj, |n| . 



However, little is known about the dynamical process 
that results in inflaton decay. We may invert this prob- 
lem, and ask the question as to how the inflaton field 
may be generated coherently from the radiation fields? 
Although the coupling between the two is small, our 
view is that Bose enhancement of the inflaton in the fi- 
nal state can lead to coherence, similar to a condensate 
effect. Incorporating this feature in the Boltzmann equa- 
tions that describe the perturbative decay of the inflaton 
is not hard. We may run the system of equations (see be- 
low) backwards in time (decreasing scale factor) to study 
how the energy densities evolve. We show that an ex- 
ponential enhancement of the inflaton density can over- 
come the Boltzmann weight, very similar to the Hagedorn 
picture. The temperature at this instant is the limiting 
temperature, beyond which the radiation energy density 
drops rapidly to zero. This may be interpreted as T max , 
the maximum or limiting temperature. Thus, we have 
an alternate physical description of the maximum tem- 
perature, and we will quantify it both theoretically and 
numerically in following. 

We consider the interaction Li nt — g<t>ff where / de- 
notes a fermion 1 . The evolution equation for the inflaton 
energy density p^ is given by 



+ 'iHp (j} = r 0( O0(l + n B {m$)) , 



(1) 



r = -wz— = a<t> m 4> . ( 2 ) 
oZir 

where H is the Hubble expansion rate and is the in- 
flaton decay rate. Note that the source term for p^ is 
enhanced by the Bose factor (I + n B (m^)) = 1/(1 — 
e _m */ T ), since the inflaton is a scalar field. Since the 
expression for is derived for zero total momentum 
(center-of-mass frame), the p = mode of the inflaton is 
being populated, p^ as defined here is the thermal energy 
density 



(3) 



(e m */ T - 1) 



The thermal factor comes out naturally when one con- 
siders fcrmions described by a thermal distribution fusing 
to form the inflaton. The analogous equation for pn, the 
energy density of radiation is 



Pr + 4:Hpn = -r^(l + n B {m,p)) 



(4) 



We evolve eqns. Q and (J2J) numerically, for which pur- 
pose it is convenient to define rescaled quantities as 



x = a/m < f > , $ = p ( pa 3 1 R = p R a 4 / . 



(5) 



1 For Standard model quarks or leptons, such an interaction cannot 
be written down since it would not be a SU(2) gauge invariant 
monomial. For the sake of exemplification, we assume that these 
fcrmions are pure singlet. 



where a is the scale factor of the Friedman-Robertson- 
Walker metric. In terms of rescaled quantities, the evo- 
lution equations take the simple form 



dx J 



dR 

dx 



VxWTr 

x 2 <S> 



VxW+R 



(l + n B (x,R)), (6) 
{l + n B (x,R)), (7) 



g^a^,Mpim°J z . These equations are simi- 



where c = J ^g^Mpim 5 / 2 . 
lar in form to those in [15|, but the rescaling is different 
since we wish to study the build-up rather than the de- 
cay of the inflaton, for which time runs backwards. Fur- 
thermore, we include the Bose enhancement factor that 
was not required in |15|. For the numerical analysis, we 
choose Tn^ = 10 13 GeV, T rch - 10 9 GeV. The traditional 
expression for the reheat temperature corresponding to 
a thermal bath dominated solely by relativistic particles 
is given by 



rch 



V 9* 



1/4 



(8) 



where g* accounts for the relativistic degrees of freedom 
in the plasma. Using this relation, and with = 200, 
we obtain a = 2.3 x 10~ 13 and c = 8.7 x 10 37 5 . We 
begin the evolution from T rc h ~ 10 9 GeV, which fixes the 
radiation energy density through the Stefan-Boltzmann 
law 
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T 4 « 2.85 x 10 36 GeV 4 



(9) 



However note that in our case there is no inflaton oscil- 
lation dominated phase and radiation is the only compo- 
nent we start with. In this respect the traditional notion 
of reheat temperature, which applies to a largest tem- 
perature of a thermal bath only dominated by relativistic 
particles, does not hold in this scenario. Nevertheless the 
numerical significance of eqn. JSJl continues to be impor- 
tant as an estimation of a temperature of a relativistic 
bath when the entire inflaton quanta has been converted 
into radiation. 

The qualitative features of the result of the numeri- 
cal analysis are independent of the choice of initial value 
of scale factor a, although x, <&, R are quantitatively de- 
pendent on it. We choose a = 10 s (x = 10 -5 ) at the 
beginning, and the results of the evolution across nearly 
4 e-foldings are displayed in Fig. ^ 

While the radiation curve R(x) (long dashed line) is re- 
maining constant, the temperature T cx a -1 , as expected 
for the radiation-dominated epoch. This is because both 
Pr and a are changing. As the inflaton energy density 
(solid line) builds up, the radiation energy density drops 
rapidly to zero. At this point, it is easy to see from 
eqn.© that <&(x) cx x 3 , which is the usual initial condi- 
tion obtained in standard evolution of the inflaton decay 
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FIG. 1: Evolution (backwards in time or with decreasing scale 
factor) of the inflaton and radiation energy densities. The 
value of the limiting temperature is at the peak T max pb m^. 



equation. Physically, our picture is quite different from 
that in |l5|. in that we can interpret T max as a limiting 
temperature, where the inflaton fields acquires a coherent 
nature, and the phase transition to an inflaton-dominated 
universe is underway. This interpretation, though novel, 
can be supported by a simple analytical estimate of T max 
which we presently derive. 

It is also interesting to observe that in our case the 
process of thermalization is actually occurring very close 
to T max in such a way that entropy density of a radiation 
dominated bath follows an adiabatic expansion law so? — 
constant, where s is the entropy density. In this respect 
we are tempted to associate T max as an actual reheat 
temperature of the universe. 

At the limiting temperature, when the inflaton energy 
density is dominant, H holds, and we may neglect 

the iHfjfj, term in eqn. . Assuming the temperature to 
be approximately constant over a small range of time, we 
may integrate the equation from time to to t (to < t in 
the absolute sense) to obtain 



In 



n^to)) ~ 1 - e -"V T 



where = p^/m^. This equation governs the rate of 
increase in the inflaton number density in a small interval 
t — to near the transition. This time interval must be 
of the order of (or slightly less than) 1/1^ < 1 / (a^m^) , 
because subsequently, the inflaton density drops off much 
more rapidly than the temperature. This means we may 
approximate as 



Now, since n^{to) includes a statistical weight factor, 
g-m^/T^ j- we assume < m ^ so that it is approximately 
a Boltzmann weight), the growth in number density of in- 
flatons begins to overwhelm the exponential suppression 
from the statistical weight when 



T sa 7710(1 



t>/T\ 



(12) 



n^t) w n^,(t )e 



i/(i- 



which is uniquely satisfied for T « O.74tti0. We re- 
gard this as a maximum or limiting temperature, reminis- 
cent of the Hagedorn temperature from hadronic physics. 
This is the main result of this work. Beyond this temper- 
ature, the radiation fields overlap strongly and populate 
1 o~ 6 the zero momentum mode of the inflaton at a rate that is 

sufficient to overcome the thermal suppression. The uni- 
verse then has only an inflaton component to its energy 
density. We see that the analytical estimate of T max or 
7h is close to the numerical value obtained from evolving 
the Boltzmann equations for inflaton decay backwards in 
time. This numerical estimate is robust and is not af- 
fected by inclusion of inverse reaction terms in the colli- 
sion kernel of the Boltzmann equations, which only be- 
come significant at temperatures exceeding the limiting 
temperature. In modulo extra assumptions on the 
gauge interactions of the fermions and their coupling to 
the gauge fields, it was argued that T max could exceed 
m<f,. We find that in the Hagedorn picture, including 
thermal masses for the fermions, as in |l6j . does not af- 
fect the bound T max w 0.74m^. 

In summary, we have shown that it is possible to inter- 
pret the maximum temperature reached in the evolution 
of inflaton decay as a Hagedorn temperature. Viewed 
backwards in time, or with decreasing size of the uni- 
verse, at this temperature, which is close to m^, coher- 
ent population of the zero momentum state of the in- 
flaton takes place. The rate of population exceeds the 
thermal suppression factor, and it is not appropriate to 
use the evolution equations any longer, since a clear dis- 
tinction between radiation and inflaton fields cannot be 
made. This is similar to the purported hadron gas to 
quark-gluon plasma (QGP) transition at finite tempera- 
ture T ~ m^, where the Hagedorn temperature signifies 
when the quark fields start to percolate through the over- 

(10) lapping hadrons. At this point, hadrons are no longer 
genuine physical degrees of freedom, and the deconfinc- 
ment phase transition is approached. We have applied 
the same idea to the early stage of inflaton decay, and 
supported the conjecture with a numerical analysis of 
the Boltzmann equations. 

We find that the maximum temperature, or the corre- 
sponding Hagedorn temperature of the post-inflationary 
universe is very close to the inflaton mass, T max fs 
O.74m0. We regard this is as an actual reheat tempera- 
ture of the universe because below the phase transition, 

(11) the subsequent evolution of radiation density obeys an 
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adiabatic expansion law with Toe a -1 . We consider this 
result as having an impact on particle-cosmology which 
we shall explore in a separate publication. 
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